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Abstract. Measurements of luminosity obtained using the ATLAS detector during early running of the Large Hadron 
Collider (LHC) at ^/l = 7 TeV are presented. The luminosity is independently determined using several detectors and 
multiple algorithms, each having different acceptances, systematic uncertainties and sensitivity to background. The 
ratios of the luminosities obtained from these methods are monitored as a function of time and of fl, the average 
number of inelastic interactions per bunch crossing. Residual time- and /i-dependence between the methods is less 
than 2% for < fl < 2.5. Absolute luminosity calibrations, performed using beam separation scans, have a common 



^ ' systematic uncertainty of ±11%, dominated by the measurement of the LHC beam currents. After calibration, the 

OO . luminosities obtained from the different methods differ by at most ±2%. The visible cross sections measured using the 

beam scans are compared to predictions obtained with the PYTHIA and PHOJET event generators and the ATLAS 
' detector simulation. 



PACS. 07.77.Ka Charged-particle beams, sources and detectors - 29.27. -a Charged-particle beams in accelerators - 
5_( ' 13.75.Cs, 13.85.-t Proton-proton interactions 
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1 Introduction and Overview 

A major goal of the ATLAS (T\ physics program for 2010 is the 
measurement of cross sections for Standard Model processes. 
Accurate determination of the luminosity is an essential ingre- 
dient of this program. This article describes the first results on 
luminosity determination, including an assessment of the sys- 
tematic uncertainties, for data taken at the LHC 12] in proton- 



proton collisions at a centre-of-mass energy — 7 TeV. It is 
organized as follows. 

The ATLAS strategy for measuring and calibrating the lu- 
minosity is outlined below and is followed in Section |2] by a 
brief description of the subdetectors used for luminosity de- 
termination. Each of these detectors is associated with one or 
more luminosity algorithms, described in Section [3] The ab- 
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solute calibration of these algorithms using beam-separation 
scans forms the subject of Section|4] The internal consistency 
of the luminosity measurements is assessed in Section |5] Fi- 
nally, the scan-based caUbrations are compared in Section |6] 
to those predicted using the PYTHIA 13] and PHOJET H 
event generators coupled to a full GEANT4 ^ simulation of 
the ATLAS detector response [6 |. Conclusions are summarized 
in Section |7] 

The luminosity of a pp collider can be expressed as 



where 7?,„g/ is the rate of inelastic collisions and (7,„g/ is the 
pp inelastic cross section. If a collider operates at a revolution 
frequency /, and n^, bunches cross at the interaction point, this 
expression can be rewritten as 

(2) 

where jj. is the average number of inelastic interactions per 
bunch crossing (BC). Thus, the instantaneous luminosity can 
be determined using any method that measures the ratio / Oi„ei- 
A fundamental ingredient of the ATLAS strategy to assess 
and control the systematic uncertainties affecting the absolute 
luminosity determination is to compare the measurements of 
several luminosity detectors, most of which use more than one 
counting technique. These multiple detectors and algorithms 
are characterized by significantly different acceptance, response 
to pile-up (multiple pp interactions within the same bunch cross- 
ing), and sensitivity to instrumental effects and to beam-induced 
backgrounds. The level of consistency across the various meth- 
ods, over the full range of single-bunch luminosities and beam 



conditions, provides valuable cross-checks as well as an esti- 
mate of the detector-related systematic uncertainties. 

Techniques for luminosity determination can be classified 
as follows: 

- Event Counting: here one determines the fraction of bunch 
crossings during which a specified detector registers an 
"event" satisfying a given selection requirement. For in- 
stance, a bunch crossing can be said to contain an "event" if 
at least one pp interaction in that crossing induces at least 
one observed hit in the detector being considered. 

- Hit Counting: here one counts the number of hits (for ex- 
ample, electronic channels or energy clusters above a spec- 
ified threshold) per bunch crossing in a given detector 

- Particle Counting: here one determines the distribution of 
the number of particles per beam crossing (or its mean) 
inferred from reconstructed quantities {e.g. tracks), from 
pulse-height distributions or from other observables that re- 
flect the instantaneous particle flux traversing the detector 
{e.g. the total ionization current drawn by a liquid-argon 
calorimeter sector). 

At present, ATLAS relies only on event-counting methods 
for the determination of the absolute luminosity. Equation |2] 
can be rewritten as: 

where £ is the efficiency for one inelastic pp collision to sat- 
isfy the event-selection criteria, and /x'" = efi is the average 
number of visible inelastic interactions per BC {i.e. the mean 
number of pp collisions per BC that pass that "event" selec- 
tion). The visible cross section (7,.(i = e<7,„e/ is the calibration 



constant that relates the measurable quantity jj."' to the lumi- provide an absolute luminosity calibration at ATLAS through 
nosity Both £ and a,,,5 depend on the pseudorapidity distri- the measurement of elastic ppscattering at small angles in the 
bution and particle composition of the collision products, and Coulomb-Nuclear Interference region. In addition, it is possi- 
are therefore different for each luminosity detector and algo- ble to normalize cross section measurements to electroweak 
rithm. processes for which precise NNLO calculations exist, for ex- 

In the limit /Xv,( ^ 1, the average number of visible inelastic ample W and Z production ifTOl . Although the cross section 
interactions per BC is given by the intuitive expression for the production of electroweak bosons in pp collisions at 

A' a/s = 7 TeV has been measured by ATLAS 1 1 1] and found to 

be in agreement with the Standard Model expectation, with ex- 

where is the number of events passing the selection criteria 

perimental and theoretical systematic uncertainties of ^ 7%, 

that are observed during a given time interval, and Nbc is the 

we choose not to use these data as a luminosity calibration, 

number of bunch crossings in that same interval. When ji in- 

since such use would preclude future comparisons with theory. 

creases, the probability that two or more pp interactions occur 

However, in the future, it will be possible to monitor the varia- 

in the same bunch crossing is no longer negligible, and /i'" 

tion of luminosity with time using W and Z production rates. 

is no longer Unearly related to the raw event count A^. Instead 

An alternative is to calibrate the counting techniques using 

/x''* must be calculated taking into account Poisson statistics, 

the absolute luminosity ^ inferred from measured accelerator 

and in some cases, instrumental or pile-up related effects (Sec- 

parameters 11121 0131 : 



tion 

Several methods can be used to determine ff,,w- At the Teva- ini^x^ 

ti-on, luminosity measurements are normalized to the total in- ^^ere m and «2 are the numbers of particles in the two collid- 
elastic pp cross section, with simulated data used to determine ^'^^'^^es and Z, and Z,, characterize the widths of the hor- 

the event- or hit-counting efficiencies Ql. Unlike the case of ^^"^^1 ^"^^ ^^^^'^^1 ^^^"^ P™^!^^- 'ypically measures 

the Tevatron, where the pp cross section was determined!] in- ^"^^ ^^^"8 ^^r Meer (vdM) scans (sometimes also called 

dependents by two experiments, the pp inelastic cross section beam-separation or luminosity scans) fll ■ The observed event 

at 7 TeV has not been measured yet. Extrapolations from lower '^'^ '^"""'^^^ ^^ile scanning the two beams across each other 

energy involve significant systematic uncertainties, as does the ^''^^ horizontal (x), then in the vertical (y) direction. This 

determination of e, which depends on the modeling of parti- measurement yields two bell-shaped curves, with the maxi- 

cle momentum distributions and multipHcity for the full pp in- "^""^ '^'^ separation, from which one extracts the values 

elastic cross section. In the future, the ALFA detector 19J will ^"^^ ^> (Section©. The luminosity at zero separation can 

' In fact, Tevatron cross sections were measured at ^= 1.8 TeV then be computed using Equation |5] and Oyis extracted from 

and extrapolated to = 1 .96 TeV. Equation |3]using the measured values of if and pi"'. 



The vdM technique allows the determination of Ovis with- rithms are added. This infrastructure enables comparison of the 

out a priori knowledge of the inelastic pp cross section or of results from different methods as a function of time. After data 

detector efficiencies. Scan results can therefore be used to test quality checks have been performed and calibrations have been 

the reliability of Monte Carlo event generators and of the AT- validated, one algorithm is chosen as the "preferred" offline al- 

LAS simulation by comparing the visible cross sections pre- gorithm for physics analysis and stored as such in the database. 

dieted by the Monte Carlo for various detectors and algorithms Luminosity information is stored as delivered luminosity. Cor- 

to those obtained from the scan data. rections for trigger prescales, DAQ deadtime and other sources 

\rrj AC jj^ J t u. • u 1 . 1 • of data loss are performed on an LB-by-LB basis when the in- 

ATLAS uses the vdM method to obtam its absolute lumi- ^ -' 

i u u f 1- • J c ca- tegrated luminosity is calculated, 

nosity calibration both for online monitoring and for offline ■' 

analysis. Online, the luminosity at the ATLAS interaction point 

(IPl) is determined approximately once per second using the 2 The ATLAS LuiTlinOSity DeteCtOrS 

counting rates from the detectors and algorithms described in 

The ATLAS detector is described in detail in Ref. |[T1. This sec- 
Sections 12] and [3] The raw event count is converted to a 

tion provides a brief description of the subsystems used for In- 
visible average number of interactions per crossing /Xjn as de- 

minosity measurements, arranged in order of increasing pseu- 

scribed in Section \3A[ and expressed as an absolute luminos- □ 

dorapidityo A summary of the relevant characteristics of these 

ity using the visible cross sections (7,,(i measured during beam- 
detectors is given in Table [T] 

separation scans. The results of all the methods are displayed 

The Inner Detector is used to measure the momentum of 

in the ATLAS control room, and the luminosity from a single 

charged particles. It consists of three subsystems: a pixel de- 

online "preferred" algorithm is transmitted to the LHC control 

tector, a silicon strip tracker (SCT) and a transition radiation 

room, providing real-time feedback for accelerator tuning. 

straw tube tracker (TRT). These detectors are located inside a 

The basic time unit for storing luminosity information for 

solenoidal magnet that provides a 2 T axial field. The tracking 

later use is the Luminosity Block (LB). The duration of a LB is — 

ATLAS uses a coordinate system where the nominal interaction 

approximately two minutes, with begin and end times set by _ 

point is at the centre of the detector The direction of beam 2 (coun- 

the ATLAS data acquisition system (DAO). All data-quality , , ■ j , t tt^ ■ x j ^ i ■ u i 

^ J V -1 J terclockwise around the LHC ring) defines the z-axis; the x-y plane is 

information, as well as the luminosity, are stored in a rela- transverse to the beam. The positive x-axis is defined as pointing to 

tional database for each LB. The luminosity tables in the offline t^e centre of the ring, and the positive j-axis upwards. Side-A of the 

database allow for storage of multiple methods for luminos- detector is on the-positive z side and side-C on the negative-z side. The 

ity determination and are versioned so that updated calibration azimuthal angle (j) is measured around the beam axis. The pseudora- 

constants can be applied. The results of all online luminosity pidity tj is defined as r] = — ln(tan6/2) where 9 is the polar angle 

methods are stored, and results from additional offline algo- from the beam axis. 



Detector 


Pseudorapidity Coverage 


# Readout Channels 


Pixel 


|r]| <2.5 


8x lO'^ 


SCT 


|ri| <2.5 


6.3 X 10* 


TRT 


|J]|<2.0 


3x 10^ 


MBTS 


2.09<|t)| <3.84 


32 


LAr: EMEC 


2.5<|t]| <3.2 


3x lO'* 


LAr: FCal 


3.1 < |t)| <4.9 


5632 


BCM 


hi =4.2 


8 


LUCID 


5.6<|t)| <6.0 


32 


ZDC 


hi > 8.3 


16 



Table 1. Summary of relevant characteristics of the detectors used 
for luminosity measurements. For the ZDC, the number of readout 
channels only includes those used by the luminosity algorithms. 



efficiency as a function of transverse momentum (pr), aver- 
aged over all pseudorapidity, rises from ^ 10% at 100 MeV to 
~ 86% for pT above a few GeV ifBll . 

For the initial running period at low instantaneous lumi- 
nosity (< 10^'^ cm^^s^^), ATLAS has been equipped with seg- 
mented scintillator counters, the Minimum Bias Trigger Scin- 
tillators (MBTS), located at z = ±365 cm from the collision 
centre. The main purpose of the MBTS is to provide a trig- 
ger on minimum collision activity during a pp bunch crossing. 
Light emitted by the scintillators is collected by wavelength- 
shifting optical fibers and guided to a photomultiplier tube 
(PMT). The MBTS signals, after being shaped and amplified, 
are fed into leading-edge discriminators and sent to the central 
trigger processor (CTP). An MBTS hit is defined as a signal 
above the discriminator threshold (50 mV). 



The precise timing (^ 1 ns) provided by the liquid argon 
(LAr) calorimeter is used to count events with collisions, there- 
fore providing a measurement of the luminosity. The LAr calor- 
imeter covers the region \rj\ < 4.9. It consists of the electro- 
magnetic (EM) for hi < 3.2, the Hadronic Endcap for 1.5 < 
|tj| < 3.2 and the Forward Calorimeter (FCal) for 3.1 < |tj| < 
4.9. The luminosity analysis is based on energy deposits in the 
Inner Wheel of the electromagnetic endcap (EMEC) and the 
first layer of the FCal. The precise timing is used to reject back- 
ground for the offline measurement of the luminosity. 

The primary purpose of the Beam Conditions Monitor 
(BCM) lfT6l is to monitor beam losses and provide fast feed- 
back to the accelerator operations team. It is an essential ingre- 
dient of the detector protection system, providing a fast acceler- 
ator abort signal in the event of large beam loss. The BCM con- 
sists of two arms of diamond sensors located at z = ±184 cm 
and r = 5.5 cm and uses programable front-end electronics 
(FPGAs) to histogram the single-sided and coincidence rates 
as a function of Bunch Crossing Identifier (BCID). These his- 
tograms are read out by the BCM monitoring software and 
made available to other online applications through the online 
network. Thus, bunch-by-bunch rates are available and are not 
subject to DAQ deadtime. The detector's value as a luminos- 
ity monitor is further enhanced by its excellent timing ( 0.7 ns) 
which allows for rejection of backgrounds from beam-halo. 

LUCID is a Cherenkov detector specifically designed for 
measuring the luminosity in ATLAS. Sixteen optically reflect- 
ing aluminum tubes filled with C4Fi() gas surround the beampipe 
on each side of the interaction point. Cerenkov photons created 
by charged particles in the gas are reflected by the tube walls 
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until they reach PMTs situated at the back end of the tubes. 
The Cherenkov light created in the gas typically produces 60- 
70 photoelectrons, while the quartz window adds another 40 
photoelectrons to the signal. After amplification, the signals 
are split three-fold and presented to a set of constant fraction 
discriminators (CFDs), charge-to-digital converters and 32-bit 
flash ADCs with 80 samplings. If the signal has a pulse height 
larger than the discriminator threshold (which is equivalent to 
15 photoelectrons) a tube is "hit." The hit-pattern produced 
by all the discriminators is sent to a custom-built electronics 
card (LUMAT) which contains FPGAs that can be programmed 
with different luminosity algorithms. LUMAT receives timing 
signals from the LHC clock used for synchronizing all detec- 
tors and counts the number of events or hits passing each lu- 
minosity algorithm for each BCID in an orbit. It also records 
the number of orbits made by the protons in the LHC during 
the counting interval. At present there are four algorithms im- 
plemented in the LUMAT firmware (see Section 13.2.31 ). The 
data from LUMAT are broadcast to the ATLAS online network 
and archived for later offline use. In addition, LUMAT provides 
triggers for the CTP and sends the hit-patterns to the DAQ. The 
LUCID electronics is decoupled from the DAQ so that it can 
provide an online luminosity determination even if no global 
ATLAS run is in progress. 

The primary purpose of the Zero-Degree Calorimeter (ZDC) 
is to detect forward neutrons and photons with |t7| > 8.3 in 
both pp and heavy-ion collisions. The ZDC consists of two 
arms located at z = ±140 m in slots in the LHC TAN (Target 
Absorber Neutral) [21, occupying space that would otherwise 
contain inert copper shielding bars. In its final configuration. 



each arm consists of calorimeter modules, one electromagnetic 
(EM) module (about 29 radiation lengths deep) followed by 
three hadronic modules (each about L14 interaction lengths 
deep). The modules are composed of tungsten with an embed- 
ded matrix of quartz rods which are coupled to photo multi- 
plier tubes and read out through CFDs. Until July 2010 only 
the three hadronic modules were installed to allow running of 
the LHCf experiment 117 1, which occupied the location where 
the EM module currently sits. Taking into account the limiting 
aperture of the beamUne, the effective ZDC acceptance for neu- 
trals corresponds to 1 GeV in pr for a 3.5 TeV neutron or pho- 
ton. Charged particles are swept out of the ZDC acceptance by 
the final-triplet quadrupoles; Monte Carlo studies have shown 
that neutral secondaries contribute a negligible amount to the 
typical ZDC energy. A hit in the ZDC is defined as an energy 
deposit above CFD threshold. The ZDC is fully efficient for 
energies above ^ 400 GeV. 



3 Luminosity Algorithms 

The time structure of the LHC beams and its consequences 
for the luminosity measurement (Section 13. Il l drive the archi- 
tecture of the online luminosity infrastructure and algorithms 
(Section |3.2| |. Some approaches to luminosity determination, 
however, are only possible offline (Section 13.3b . In all cases, 
dealing properly with pile-up dependent effects (Section 13.4b 
is essential to ensure the precision of the luminosity measure- 
ments. 



3.1 Bunch Patterns and Luminosity Bacl<grounds 

The LHC beam is subdivided into 35640 RF-buckets of which 
nominally every tenth can contain a bunch. Subtracting abort 
and injection gaps, up to 2808 of these 3564 "slots", which 
are 25 ns long, can be filled with beam. Each of these possible 
crossings is labeled by an integer BCID which is stored as part 
of the ATLAS event record. 

Figure[T]displays the event rate per BC, as measured by two 
LUCID algorithms, as a function of BCID and time-averaged 
over a run that lasted about 15 hours. For this run, 35 bunch 
pairs collided in both ATLAS and CMS. These are called "col- 
liding" (or "paired") BCIDs. Bunches that do not colUde at IPl 
are labelled "unpaired." Unpaired bunches that undergo no col- 
lisions in any of the IPs are called "isolated." The structures 
observed in this figure are visible in the bunch-by-bunch lu- 
minosity distributions of all the detectors discussed in this pa- 
per, although with magnitudes affected by different instrumen- 
tal characteristics and background sensitivities. Comparisons 
of the event rates in colliding, unpaired, isolated and empty 
bunch crossings for different event-selection criteria provide 
information about the origin of the luminosity backgrounds, 
as well as quantitative estimates of the signal purity for each of 
these detectors and algorithms. 

Requiring at least one hit on at least one side (this is re- 
ferred to as an Event_OR algorithm below) reveals a complex 
time structure (Fig.[T^). The colliding bunches are clearly dis- 
tinguished, with a rate of about four orders of magnitude above 
background. They are followed by a long tail where the rate 
builds up when the paired BCID's follow each other in close 
succession, but decays slowly when no collisions occur for a 
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sufficiently long time. This "afterglow" (which is also appar- 
ent when analyzing the luminosity response of Event_OR al- 
gorithms using the BCM or MBTS) is dominated by slowly- 
decaying, low-energy radiation produced by pp collision prod- 
ucts that hit forward ATLAS components and scatter around 
the experimental cavern for tens of microseconds. BCID's from 
unpaired and isolated bunches appear as small spikes above the 
afterglow background. These spikes are the result of beam-gas 
and beam-halo interactions; in some cases, they may also con- 
tain a very small fraction of pp collisions between an unpaired 
bunch in one beam and a satellite- or debunched- proton com- 

y 

ponent in the opposing beamO 

For the Event_AND algorithm (Fig.[TJ)), the coincidence re- 
quirement between the A- and C-sides suppresses the afterglow 
signal by an additional four orders of magnitude, clearly show- 
ing that this luminosity background is caused by random sig- 
nals uncorrected between the two sides. Unpaired-bunch rates 
for LUCID_EventJ\ND lie 4-5 orders of magnitude lower than 
pp collisions between paired bunches. 

This figure illustrates several important points. First, be- 
cause only a fraction of the BCID's are filled, an algorithm 
that selects on colliding BCID's is significantly cleaner than 
one that is BCID-blind. Second, and provided only colliding 
BCID's are used, the background is small (LUCID) to moder- 

^ In proton storage rings, a small fraction of the injected (or stored) 
beam may fail to be captured into (or may slowly diffuse out of) 
the intended RF bucket, generating a barely detectable unbunched 
beam component and/or coalescing into very low-intensity "satellite" 
bunches that are separated from a nominal bunch by up to a few tens 
of buckets. 
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ate (MBTS) for Event.OR algorithms, and negligible for 
Event^ND. In the Event.OR case, the background contains 
contributions both from afterglow and from beam-gas and beam- 
halo interactions: its level thus depends crucially on the time 
separation between colliding bunches. 

3.2 Online Algorithms 

3.2.1 Online Luminosity Infrastructure 

Online luminosity monitoring and archiving can be made avail- 
able even when only the core ATLAS DAQ infrastructure is 
active; this makes it possible to provide luminosity informa- 
tion for machine tuning independently of the "busy" state of 
the DAQ system and of the hardware status of most subdetec- 
tors (except for the CTP and for one or more of the luminosity 
detectors). In addition, since the online luminosity data are col- 
lected in the front-end electronics of each detector (or at the 
CTP input), there is no need for prescaling, even at the highest 
luminosities. 

The calculation and publication of instantaneous luminosi- 
ties is performed by an application suite called the Online Lu- 
minosity Calculator (OLC). The task of the OLC is to retrieve 
the raw luminosity information (event or hit counts, number of 
colliding bunches rtj,, and number of LHC orbits in the time in- 
terval considered) from the online network and to use these data 
to determine jj, and hence the measured luminosity. For each 
luminosity algorithm, the OLC outputs the instantaneous lumi- 
nosity, averaged over all colliding BClDs, at about 1 Hz. These 
values are displayed on online monitors, stored in the ATLAS 
online-monitoring archive and shipped to the LHC control room 



to assist in collision optimization at IPl. In addition, the OLC 
calculates the luminosity averaged over the current luminos- 
ity block (in all cases the luminosity averaged over all collid- 
ing BClDs, and when available the bunch-by-bunch luminosity 
vector) and stores these in the ATLAS conditions database. 

Most methods provide an LB-averaged luminosity measured 
from colliding bunches only, but for different detectors the re- 
quirement is imposed at different stages of the analysis. The 
BCM readout driver and the LUCID LUMAT module provide 
bunch-by-bunch raw luminosity information for each LB, as 
well as the luminosity per LB summed over all colliding BCID's. 
For these two detectors, the OLC calculates the total (i.e. bunch- 
integrated) luminosity using an extension of Equation |3] that 
remains valid even when each bunch pair produces a different 
luminosity (reflecting a different value of jj.) because of differ- 
ent bunch currents and/or emittances: 

E (6) 

ieBCID 

where the sum is performed over the colliding BCID's. This 
makes it possible to properly apply the pile-up correction bunch- 
by-bunch (Section [34l l. 

For detectors where bunch-by-bunch luminosity is unavail- 
able online. Equation [3] is used, with jU™ computed using the 
known number of paired BCID's and the raw luminosity infor- 
mation averaged over either the colliding BCID's (this is the 
case for the MBTS) or all BCID's (the front-end luminosity in- 
frastructure of the ZDC provides no bunch-by-bunch capability 
at this time). 

For the MBTS, which lacks appropriate FPGA capabili- 
ties in the front end, the selection of colliding bunches is done 
through the trigger system. The BCID's that correspond to col- 
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(a) (b) 
Fig. 1. Bunch-by-bunch event rate per bunch crossing in ATLAS run 162882, as recorded by a LUCID algorithm that requires (a) at least one 
hit on either LUCID side (Event.OR), or (b) at least one hit on both LUCID sides (Event_AND) within the same BCID. 



liding bunches are identified and grouped in a list called the 
"physics bunch group," which is used to gate the physics trig- 
gers. A second set of triggers using unpaired bunches is used 
offline to estimate beam backgrounds. The MBTS counters pro- 
vide trigger signals to the CTP, which then uses bunch-group 
information to create separate triggers for physics and for un- 
paired bunch groups. The CTP scalers count the number of 
events that fire each trigger, as well as the number of LHC or- 
bits (needed to compute the rate per bunch crossing). Every 
10 s these scalers are read out and pubhshed to the online net- 
work. Three values are stored for each trigger type: trigger be- 
fore prescale (TBP), trigger after prescale and trigger after veto 
(TAV). The TBP counts are calculated directly using inputs to 
the CTP and are therefore free from any dead time or veto (ex- 
cept when the DAQ is paused), while the TAV corresponds to 
the rate of accepted events for which a trigger fired. To maxi- 
mize the statistical power of the measurement and remain un- 



affected by prescale changes, online luminosity measurements 
by the MBTS algorithms use the TBP rates. 

3.2.2 BCM Algorithms 

Out of the four sensors on each BCM side, only two are cur- 
rently used for online luminosity determination. Three onUne 
algorithms, implemented in the firmware of the BCM readout 
driver, report results: 

- BCM_Event_OR counts the number of events per BC in which 
at least one hit above threshold occurs on either the A-side, 
the C-side or both, within a 12.5 ns window centred on the 
arrival time of particles originating at IPl; 

- BCM_Event J^N D counts the number of events per BC where 
at least one hit above threshold is observed, within a 12.5 ns- 
wide coincidence window, both on the A- and the C-side. 
Because the geometric coverage of the BCM is quite small. 
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the event rate reported by this algorithm during the beam- 
separation scans was too low to perform a reliable calibra- 
tion. Therefore this algorithm will not be considered further 
in this paper; 

- BCM_Event_XORC counts the number of events per BC 
where at least one hit above threshold is observed on the 
C-side, with none observed on the A-side within the same 
12.5 ns-wide window. Because converting the event-counting 
probabiUty measured by this method into an instantaneous 
luminosity involves more complex combinatorics than for 
the simpler Event.OR and Event_AND cases, fully exploit- 
ing this algorithm requires more extensive studies. These 
he beyond the scope of the present paper. 

3.2.3 LUCID Algorithms 

Four algorithms are currently implemented in the LUMAT card: 

- LLICID_Zero_OR counts the number of events per BC where 
at least one of the two detector sides reports no hits within 
one BCID, or where neither side contains any hit in one 
BCID; 

- LUCID_Zero_AND counts the number of events per BC where 
no hit is found within one BCID on either detector side; 

- LUCID.HIt.OR reports the mean number of hits per BC. In 
this algorithm, hits are counted for any event where there is 
at least one hit in any one of the 16 tubes in either detector 
side in one BCID; 

- LUCID_HltJ\ND reports the mean number of hits per BC, 
with the additional requirement that the event contain at 
least one hit on each of the two detector sides in one BCID. 



The LUCID event-counting algorithms simply subtract the 
number of empty events reported by the zero-counting algo- 
rithms above from the total number of bunch crossings: 

- LUCID_Event_AND reports the number of events with at least 
one hit on each detector side (A^LUCID_Event_AND ~ ^bc — 

^LUCID_Zero.OR); 

- LUCID.Event.OR reports the number of events for which 

the sum of the hits on both detector sides is at least one 
(^LUCID.Event.OR = ^bc - A^LUCID^eroJVND)- 

Converting measured hit-counting probabiUties into instan- 
taneous luminosity does not lend itself to analytic models of the 
type used for event counting and requires detailed Monte Carlo 
modeling that depends on the knowledge of both the detector 
response and the particle spectrum in pp colUsions. This mod- 
eUng introduces additional systematic uncertainties and to be 
used reUably requires more extensive studies that he beyond 
the scope of the present paper. 

3.2.4 MBTS Algorithms 

Raw onUne luminosity information is supplied by the following 
two CTP scalers: 

- MBTS_Event_OR counts the number of events per BC where 
at least one hit above threshold is observed on either the A- 
side or the C-side, or both; 

- MBTS-Event J^ND counts the number of events per BC where 
at least one hit above threshold is observed both on the A- 
and the C-side. 



3.2.5 ZDC Algorithms 

Online luminosity information is supplied by dedicated ZDC 
scalers that count pulses produced by constant-fraction discrim- 
inators connected to the analog sum of ZDC photomultipher 
signals on each side separately: 

- ZDC_A reports the event rate where at least one hit above 
threshold is observed on the A-side, irrespective of whether 
a hit is simultaneously observed on the C-side; 

- ZDC-C reports the event rate where at least one hit above 
threshold is observed on the C-side, irrespective of whether 
a hit is simultaneously observed on the A-side; 

- ZDC_Event_AND reports the event rate where at least one 
hit above threshold is observed in coincidence on the A- 
and C-sides. This algorithm is still under study and is not 
considered further in this paper. 

The data described here were taken before the ZDC electronic 
gains and timings were fully equalized. Hence the correspond- 
ing visible cross sections for the A- and C-side differ by a few 
per cent. 

3.3 Offline Algorithms 

Some luminosity algorithms require detailed information that 
is not easily accessible online. These algorithms use data col- 
lected with a minimum bias trigger (e.g. one of the MBTS 
triggers) and typically include tighter requirements to further 
reduce backgrounds. Because such analyses can only be per- 
formed on events that are recorded by the DAQ system, they are 
statistically less powerful than the online algorithms. However, 
since the MBTS rates per BCID are not available online, offline 
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algorithms are important for these detectors for runs where the 
currents are very different from one bunch to the next. In ad- 
dition, these methods use event selection criteria that are very 
similar to final physics analyses. 

Verification that the luminosities obtained from the offline 
methods agree well with those obtained from the online tech- 
niques through the full range of relevant jj. provides an impor- 
tant cross-check of systematic uncertainties. As with the online 
measurements, the LB-averaged instantaneous luminosities are 
stored in the ATLAS conditions database. 

3.3.1 MBTS Timing Algorithm 

The background rate for events passing the MBTS_Event_AND 
trigger is a factor of about 1000 below the signal. As a result, 
online luminosity measurements from that trigger can be reli- 
ably calculated without performing a background subtraction. 
However, the signal-to-background ratio is reduced when the 
two beams are displaced relative to each other (since the signal 
decreases but the beam-induced backgrounds remain constant). 
At the largest beam separations used during the vdM scans, the 
background rate approaches 10% of the signal. While these 
backgrounds are included in the fit model used to determine 
the online MBTS luminosity calibration (see Section 14.31 ). it 
is useful to cross-check these calibrations by reanalysing the 
data with a tighter offline selection. The offline time resolu- 
tion of the MBTS is ^ 3 ns and the distance between the A- 
and C-sides corresponds to a time difference of 23 ns for par- 
ticles moving at the speed of light. Imposing a requirement 
that the difference in time measured for signals from the two 
sides be less than 10 ns reduces the background rate in the 
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MBTS.EventJ\ND triggered events to a negligible level (< lO"'*) 
even at the largest beam displacements used in the scans, while 
maintaining good signal efficiency. This algorithm is called 
MBTS-Timing. In those instances where different bunches have 
substantially different luminosities, MBTS.Timing can be used 
to properly account for the pile-up dependent corrections. 

3.3.2 Liquid Argon Algoritlnm 

The timing cut used in MBTS.TIming is only applicable to co- 
incidence triggers, where hits are seen both on the A- and C- 
sides. It is possible to cross-check the online calibration of 
the single-sided MBTS_Event_OR trigger, where the signal-to- 
background ratios are lower, by imposing timing requirements 
on a different detector. The LAr_Timing algorithm uses the liq- 
uid argon endcap calorimeters for this purpose. Events are re- 
quired to pass the MBTS.Event.OR trigger and to have signif- 
icant in-time energy deposits in both EM calorimeter endcaps. 
The analysis considers the energy deposits in the EMEC In- 
ner Wheels and the first layer of the FCal, corresponding to the 
pseudorapidity range 2.5 < | Tj | < 4.9. Cells are required to have 
an energy 5c7 above the noise level and to have E > 250 MeV 
in the EMEC or £ > 1200 MeV in the FCal. Two cells are re- 
quired to pass the selection on each of the A- and C-side. The 
time on the A-side (C-side) is then defined as the average time 
of all the cells on the A-side (C-side) that pass the above re- 
quirements. The times obtained from the A-side and C-side are 
then required to agree to better than ±5 ns (the distance be- 
tween the A- and C-sides corresponds to a time difference of 
30 ns for particles moving at the speed of light). 



3.3.3 Track-Based Algoritlnms 

Luminosity measurements have also been performed offline 
by counting the rate of events with one or more reconstructed 
tracks in the MBTS.Event.OR sample. Here, rather than impos- 
ing a timing cut, the sample is selected by requiring that one 
or more charged particle tracks be reconstructed in the inner 
detector. Two variants of this analysis have been implemented 
that differ only in the details of the track selection. 

The first method, referred to here as primary-vertex event 
counting (PrimVtx) has larger acceptance. The track selection 
and vertex reconstruction requirements are identical to those 
used for the study of charged particle multiplicities at ^ = 
7 TeV fTS). Here, a reconstructed primary vertex is required 
that is formed from at least two tracks, each with pr > 100 MeV. 
Furthermore, the tracks are required to fulfill the following qual- 
ity requirements: transverse impact parameter \do\ < 4 mm with 
respect to the luminous centroid, errors on the transverse and 
longitudinal impact parameters G{do) < 5 mm and g{zo) < 
10 mm, at least 4 hits in the SCT, and at least 6 hits in Pixel 
and SCT. 

The second analysis, referred to here as charged-particle 
event counting (ChPart), is designed to allow the comparison 
of results from ALICE, ATLAS and CMS. It therefore uses 
fiducial and pj requirements that are accessible to all three ex- 
periments. The method counts the rate of events that have at 
least one track with transverse momentum pr > 0.5 GeV and 
pseudorapidity |t7| < 0.8. The track selection and acceptance 
corrections are identical (with the exception of the |rj| < 0.8 
requirement) to those in Ref. ifTSll . The main criteria are an 
IVIBTS_Event_OR trigger, a reconstructed primary vertex with 



at least three tracks with pr > 150 MeV, and at least one track 
with pt > 500 MeV, \ri \ < 0.8 and at least 6 SCT hits and one 
Pixel hit. Data are corrected for the trigger efficiency, the effi- 
ciency of the vertex requirement and the tracking efficiency, all 
of which depend on pj and Tj . 

3.4 Converting Counting Rates to Absolute 
Luminosity 

The value of jj.™ used to determine the bunch luminosity in 
BCID / is obtained from the raw number of counts A^, and the 
number of bunch crossings Nbc, using an algorithm-dependent 
expression and assuming that: 

- the number of pp-interactions occuring in any bunch cross- 
ing obeys a Poisson distribution. This assumption drives the 
combinatorial formalism presented in Sections 13.4.11 and 
[3A2l below: 

- the efficiency to detect a single inelastic pp interaction is 
constant, in the sense that it does not change when several 
interactions occur in the same bunch crossing. This is tanta- 
mount to assuming that the efficiency e„ for detecting one 
event associated with n interactions occuring in the same 
crossing is given by 

e„ = l-(l-eir (7) 

where ei is the detection efficiency corresponding to a sin- 
gle inelastic interaction in a bunch crossing (the same def- 
inition applies to the efficiencies e*^^, e^, e'^ and e^^^ 
defined below). This assumption will be validated in Sec- 
tion [3A3l 
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The bunch luminosity is then given directly and without addi- 
tional assumptions by 

= (8) 

using the value of (7,.,v measured during beam-separation scans 
for the algorithm considered. However, providing a value for 
ji = /i''"/e = fi"'^ (Jinei / (yvis requires an assumption on the as 
yet unmeasured total inelastic cross section at 

3.4.1 Inclusive-OR Algorithms 

In the Event-OR case, the logic is straightforward. Since the 
Poisson probability for observing zero events in a given bunch 
crossing is Po(m"*) = e^^'" — e^^^^"" , the probability of ob- 
serving at least one event is 

-TEvenLOR l,M ) — JJ^ 

= l-Po(Ai™) 

Here the raw event count Nqr is the number of bunch cross- 
ings, during a given time, in which at least one pp interaction 
satisfies the event-selection criteria of the OR algorithm under 
consideration, and Nbc is the total number of bunch crossings 
during the same interval. Equation|9]reduces to the intuitive re- 
sult PEvent.oR(M'") ~ M"' whcu << 1. SoMug for /i™ in 
terms of the event-counting rate yields: 

M-^ = -ln(l-^) (10) 

3.4.2 Coincidence Algoritlnms 

For the Event_AND case, the relationship between /i^'" and is 

more complicated. Instead of depending on a single efficiency, 
4 ATLAS uses the PYTHIA value of 71.5 mb. 
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the event-counting probability must be written in terms of e"^. The value of jUvf j is then obtained by solving Equation fT2lnu- 

and e'^^^, the efficiencies for observing an event with, re- merically using the values of a^^ and <7^^^ extracted from 

spectively, at least one hit on the A-side, at least one hit on beam separation scans. The validity of this technique will be 

the C-side and at least one hit on both sides simultaneously, quantified in Section |5] 

These efficiencies are related to the Event.OR efficiency by If the efficiency is high and g'^^^ « e'* « as is the case 

gOR ^ £A_^_^c _ £AND_ Jqj. MBTS.Event-AND, Equation[TT]can be approximated by 



M"-^-«-ln(l-^] (13) 



The probability PEvenij\ND(M) of there being at least one hit 
on both sides is one minus the probability f>Zero_OR there 
being no hit on at least one side. The latter, in turn, equals the ^he ^-dependence of the probabihty function 



Event_AND 



probability that there be no hit on at least side A (Pqa = e-^^'" ), controlled by the relative magnitudes of e^ and e^^^ (or of 
plus the probability that there be no hit on at least side C (Pqc = *^ coiTesponding measured visible cross sections). This is in 
e-''''), minus the probabihty that there be no hit on either side '^""^''^^^ Event.OR case, where the efficiency Eor factors 



P — '^AND 

-< EventJ^ND (H- 1 



out of Equation [TOl 



3.4.3 Pile-up-related Instrumental Effects 



(M) ^^^^ 
I _ (^g-A"^'* _|_ e^M'^'^ _ e^^^^") The /i-dependence of the probability functions PEvent.oR and 

\ _ (e-M*:" J- g-p^^^ ^ ^'Evenij\ND IS displayed in Fig. |2] All algorithms MfMrafe at high /x, 

reflecting the fact that as the pile-up increases, the probabihty 
of observing at least one event per bunch crossing approaches 
one. Any event-counting luminosity algorithm will therefore 
lose precision, and ultimately become unusable, as the LHC lu- 
minosity per bunch increases far beyond present levels. The tol- 
erable pile-up level is detector- and algorithm-dependent: the 
higher the efficiency > e^^^^ > ef«„^ > etucio)^ the 

earlier the onset of this saturation. 

The accuracy of the event-counting formalism can be ver- 
ified using simulated data. Figure |2] (bottom) shows that the 
parameterizations of Sections [3 .4. II and 13 .4.21 deviate from the 



This equation cannot be inverted analytically. The most appro- 
priate functional form depends on the values of e^, and 

^AND 

For cases such as LUCID.EventJ\ND and BCM_Event_AND, 
the above equation can be simplified using the fact that e'^^^ < < 
E^'^, and assuming that sa e'~. The efficiencies e^'^^ and 
e''* ai-e definedhy, respectively, e'^'^" = (7^^"/(7i„ei and £^'^ = 
a^if /Oineh the average number of visible inelastic interactions 
per BC is computed as /i'" = e^^^pL. Equation [TTj then be- 
comes 



— 1 — 2e ''('^ ^l^ + e ^^^^ full simulation by ±2% at most: possible instrumental effects 

= 1 —2e ' " 1" + e ' " I not accounted for by the combinatorial formalism are predicted 




10"^ 10'^ 1 10 



Fig. 2. Fraction of bunch crossings containing a detected event for 
LUCID and MBTS algoritlims as a function of fl, the true average 
number of inelastic pp interactions per BC. The plotted points are the 
result of a Monte Carlo study performed using the PYTHIA event 
generator together with a GEANT4 simulation of the ATLAS detec- 
tor response. The curves reflect the combinatorial formalism of Sec- 
tions [3AT] and [3A2l using as input only the visible cross sections 
extracted from that same simulation. The bottom inset shows the dif- 
ference between the full simulation and the parameterization. 

to have negligible impact for the bunch luminosities achieved 
in the 2010 LHC run (0 < jU < 5). 

It should be stressed, however, that the agreement between 
the Poisson formalism and the full simulation depends criti- 
cally on the validity of the assumption, summarized by Equa- 
tion |7] that the efficiency for detecting an inelastic pp inter- 
action is independent of the number of interactions that occur 
in each crossing. This requires, for instance, that the thresh- 
old for registering a hit in a phototube (nominally 15 photo- 
electrons for LUCID) be low enough compared to the aver- 
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age single-particle response. This condition is satisfied by the 
simulation shown in Fig.|2] Repeating this simulation with the 
LUCID threshold raised to 50 photoelectrons yields systematic 
discrepancies as large as 7% between the computed and sim- 
ulated probability functions for the LUCID Event_AND algo- 
rithm. When the threshold is too high, a particle from a single 
pp interaction occasionally fails to fire the discriminator. How- 
ever, if two such particles from different pp interactions in the 
same bunch crossing traverse the same tube, they may produce 
enough fight to register a hit. This effect is called migration. 

4 Absolute Calibration Using 
Beam-Separation Scans 

The primary calibration of all luminosity algorithms is derived 
from data collected during van der Meer scans. The principle 
(Section 14.1b is to measure simultaneously the collision rate 
at zero beam separation and the corresponding absolute lumi- 
nosity inferred from the charge of the colliding proton bunches 
and from the horizontal and vertical convolved beam sizes lfT3l . 
Three sets of beam scans have been carried out in ATLAS, as 
detailed in Section l4~2l These were performed in both the hor- 
izontal and the vertical directions in order to reconstruct the 
transverse convolved beam profile. During each scan, the colli- 
sion rates measured by the luminosity detectors were recorded 
while the beams were moved stepwise with respect to each 
other in the transverse plane. 
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4.1 Absolute Luminosity from Beam Parameters by the equation: 

1 jR,{5)d5 

In terms of colliding-beam parameters, the luminosity ^ is de- — , — — TTTT^ — 

v27r t<x[U) 

fined (for beams that collide with zero crossing angle) as 

In the case where the luminosity curve Rx{5) is Gaussian, Ex 
^ = ni, frn\n2 J pi{x,y)p2{x,y)dxdy (14) coincides with the standard deviation of that distribution. By 

using the last two equations, Equation[T5]can be rewritten as 

where rtj, is the number of colliding bunches, /, is the ma- 

chine revolution frequency (11 245. 5 Hz for LHC), n\{2) is the = — — (18) 

which is a general formula to extract luminosity from machine 
parameters by performing a beam separation scan. EquationfTSl 
is quite general; Zx and Zy only depend on the area under the 
luminosity curve. 

4.2 Luminosity-Scan Data Sets 



number of particles per bunch in beam 1 (2) and P\[2){x,y) is 
the normalized particle density in the transverse (x-y) plane of 
beam 1 (2) at the IP. Under the general assumption that there 
is no correlation between x and y, i.e. that transverse coupling 
is negligiblq^ at the IP, the particle densities can be factorized 
(p{x^y) ~ p{x)p{y)) and Equation fl4lre written as 



^ = nbfrnxn2 £lx{pi{x),P2{x)) £ly{pi{y),P2{y)) (15) 
where 



Three van der Meer scans have been performed at the ATLAS 
interaction point (Table |2]i. The procedure 11211191 ran as fol- 

^xiPi.Pi) = j P\{x)p2{x)dx j^^g ^jjgj. (.g„tj.i„g the beams on each other at the IP in both 

is the beam overlap integral in the x direction (with an analo- the horizontal and the vertical plane using mini-scans, a full 

gous definition in the y direction). In the method proposed by luminosity-calibration scan was carried out in the horizontal 

van der Meer fIT] the overlap integral (for example in the x plane, spanning a range of ±6c7,, in horizontal beam-separation 

direction) can be calculated as: (where is the nominal transverse size of either beam at the 

„ IP). A full luminosity-calibration scan was then carried out in 

^x{PuP2)= r/(J,. (16) 

J Kx[0)ao jjjg vertical plane, again spanning a range of ±6(7h in relative 

where Rx{5) is the luminosity (or equivalently /x''*) - at this beam separation. 

stage in arbitrary units - measured during a horizontal scan The mini-scans used to first centre the beams on each other 

at the time the two beams are separated by the distance 5 and in the transverse plane were done by activating closed orbit 

5 = represents the case of zero beam separation. Ex is defined bump^ around the IP that vary the IP positions of both beams 



Combining the observed vertical to horizontal emittance ratios 6 A closed orbit bump is a local distortion of the beam orbit that is 
with the measured LHC lattice functions indicates that the luminosity implemented using pairs of steering dipoles located on either side of 
loss caused by the residual tilt of the two beams was less than 0.25%. the affected region. In this particular case, these bumps are tuned to 



by iloi, in opposite directions, either horizontally or verti- 
cally. The relative positions of the two beams were then ad- 
justed, in each plane, to achieve (at that time) optimum trans- 
verse overlap. 

The full horizontal and vertical scans followed an identi- 
cal procedure, where the same orbit bumps were used to dis- 
place the two beams in opposite directions by ±3a/„ result- 
ing in a total variation of ±6(7/, in relative displacement at the 
IP. In Scan 1, the horizontal scan started at zero nominal sep- 
aration, moved to the maximum separation in the negative di- 
rection, stepped back to zero and on to the maximum positive 
separation, and finally returned to the original settings of the 
closed-orbit bumps (zero nominal separation). The same pro- 
cedure was followed for the vertical scan. In Scan II and III, af- 
ter collision optimization with the transverse mini-scans, a full 
horizontal scan was taken from negative to positive nominal 
separation, followed by a hysteresis cycle where the horizontal 
nominal separation was run to —6(7;,, then then +6(7f„ and 
finally followed by a full horizontal scan in the opposite direc- 
tion to check for potential hysteresis effects. The same proce- 
dure was then repeated in the vertical direction. 

For each scan, at each of 27 steps in relative displacement, 
the beams were left in a quiescent state for ^ 30 seconds. Dur- 
ing this time the (relative) luminosities measured by all ac- 
tive luminosity monitors were recorded as a function of time 
in a dedicated online-data stream, together with the value of 
the nominal separation, the beam currents and other relevant 
accelerator parameters transmitted to ATLAS by the accelera- 

translate either beam parallel to itself at the IP, in either the horizontal 
or the vertical direction. 
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tor control system. In addition, the full data acquisition system 
was operational throughout the scan, using the standard trigger 
menu, and triggered events were recorded as part of the normal 
data collection. 

4.3 Parametrization and Analysis of the Beam Scan 
Data 

Data from all three scans have been analyzed both from the 
dedicated online-data stream and from the standard ATLAS 
data stream. Analyses using the standard data stream suffer 
from reduced statistical precision relative to the dedicated stream, 
but allow for important cross-checks both of the background 
rates and of the size and position of the luminous region. In 
addition, because this stream contains full events, these data 
can be used to measure the visible cross section corresponding 
to standard analysis selections that require, for example, tim- 
ing cuts in the MBTS or the liquid argon Calorimeter or the 
presence of a reconstructed primary vertex. Measurements per- 
formed using these two streams provide a consistent interpre- 
tation of the data within the relevant statistical and systematic 
uncertainties. 

In all cases, the analyses fit the relative variation of the 
bunch luminosity as a function of the beam separation to ex- 
tract Ex and Ey (Equation [TtI i. These results are then com- 
bined with the measured bunch currents to determine the ab- 
solute luminosity using Equation [18] Although the pile-up ef- 
fects remained relatively weak during these scans, the raw rates 
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vdM Scan I 


vdM Scan II, III 




(April 26, 2010) 


(May 9, 2010) 


LHC Fill Number 


1059 


1089 


Scan Directions 


1 horizontal scan 


2 horizontal scans 




followed by 1 vertical scan 


followed by 2 vertical scans 


Total Scan Steps per Plane 


27 


54 (27+27) 




{±6ah) 


(±6(Ti) 


Scan Duration per Step 


30 sec 


30 sec 


Number of bunches colliding in ATLAS 


1 


1 


Total number of bunches per beam 


2 


2 


Number of protons per bunch 


~0.1- 10" 


~0.2- lOi' 


j8*(m) 


~2 


-2 


tTi (flm) [assuming nominal emittances] 


~45 


~45 


Crossing angle (/xrad) 








Typical luminosity /bunch (/ib^' /.s) 


4.5-10-3 


1.8-10-2 


fl (interactions/crossing) 


0.03 


0.11 



Table 2. Summary of the main characteristics of the three beam scans performed at the ATLAS interaction point. The values of luminosity /bunch 
and fl are given for zero beam separation. 



(^Event OR' ^Event_AND'-- ) converted|_|into a mean num- Here {nin2)meas is the product of the numbers of protons in the 

ber of interactions per crossing /i'" as described in Section[331 two colliding bunches during the measurement, {n\n2)MAX is 

In addition, to remove sensitivity to the slow decay of the beam its maximum value during the scans, and Rmeas is the value of 

currents over the duration of the scan, the data are analyzed /i'" at the current scan point. 



as specific rates, obtained by dividing the measured average 
interaction rate per BC by the product of the bunch currents 
measured at that scan point: 



R 



(ni«2) 



MAX 



sp 



{n\n2)meas 



R„ 



(19) 



Beam cuiTents are measured using two complementary LHC 
systems ll20l . The fast bunch-current transformers (FBCT) are 
AC-coupled, high-bandwidth devices which use gated electron- 
ics to perform continuous measurements of individual bunch 
charges for each beam. The Direct-Current Current Transform- 



^ For the coincidence algorithms, the procedure is iterative because ers (DCCT) measure the total circulating intensity in each of 
it requires the a priori knowledge of CT,,,-.,. Monte Carlo estimates were the two beams iiTespective of their underlying time structure. 



used as the starting point. 



The DCCT's have intrinsically better accuracy, but require av- 
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eraging over hundreds of seconds to achieve the needed pre- 
cision. The relative (bunch-to-bunch) currents are based on the 
FBCT measurement. The absolute scale of the bunch intensities 
n 1 and «2 is determined by rescaling the total circulating charge 
measured by the FBCTs to the more accurate DCCT measure- 
ments. Detailed discussions of the performance and caUbration 
of these systems are presented in Ref. ET\ . 

Fits to the relative luminosity require a choice of parametriza- 
tion of the shape of the scan curve. For all detectors and algo- 
rithms, fits using a single Gaussian or a single Gaussian with 
a flat background yield unacceptable distributions. In all 
cases, fits to a double Gaussian (with a common mean) plus 
a flat background result in a per degree of freedom close to 
one. In general, the background rates are consistent with zero 
for algorithms requiring a coincidence between sides, while 
small but statistically significant backgrounds are observed for 
algorithms requiring only a single side. These backgrounds are 
reduced to less than 0.3% of the luminosity at zero beam sepa- 
ration by using data from the paired bunches only. Offline anal- 
yses that require timing or a primary vertex, in addition to being 
restricted to paired bunches, have very low background. The 
residual background is subtracted using the rate measured in 
unpaired bunches; no background term is therefore needed in 
the fit function for the offline case. Examples of such fits are 
shown in Fig. [3] 

For these fits the specific rate is described by a double Gaus- 



Here a, and a, are the widths of first and second Gaussians 
respectively, /, is the fraction of the rate in the first Gaussian 
and xq is introduced to allow for the possibility that the beams 
are not perfectly centred at the time of the scan. The value of 
Ex in Equation [18] is calculated as 



1 



f. I l-fi 

CJ, Gi 



(21) 



sian: 



Rx{5)d5 



1% 



(''--^o)~ 



(20) 



4.4 Fit Results 

Summaries of the relevant fit parameters for the three scans are 
presented in Tables |7] through |9] in Appendix A. Because the 
emittance during Scan 1 was different from that during Scans II 
and III, the values of and are not expected to be the 
same for the first and the later scans. Furthermore, because the 
beam currents were lower in Scan I, the peak luminosities for 
this scan are lower than for the later scans. These tables, as 
well as Fig. |4] show that the mean position and Z for a given 
scan are consistent within statistical uncertainties amongst all 
algorithms. These data also indicate several potential sources 
of systematic uncertainty. First, the fitted position of the peak 
luminosity deviates from zero by as much as 7 /im, indicating 
that the beams may not have been properly centred before the 
start of the scan. Second, in scans II and III, the peak luminosi- 
ties for the horizontal and vertical scans, as measured with a 
single algorithm, show a systematic difference of as much as 
5% (with a lower rate observed in the vertical scan for all al- 
gorithms). This systematic dependence may indicate a level of 
irreproducibility in the scan setup. The effect of these system- 
atic uncertainties on the luminosity calibration is discussed in 
Section E 
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Fig. 3. Results of fits to the second luminosity scan in the x (left) and y (right) direction for the (a) LUCID_Event_OR, (b) MBTS.TIming, and 
(c) ChPart algorithms. The panels at the bottom of each graph show the difference of the measured rates from the value predicted by the fit, 
normalized to the statistical uncertainty on the data (<7). 



Calibration of the absolute luminosity from the beam scans 
uses the following expression for Ovis'. 

= ^MAX = ^ n,Mn,n2)MAX ^^^^ 
where r'^^^ and are, respectively, the value of Rs,, and 

the absolute luminosity (inferred from the measured machine 
parameters) when the beams collide exactly head-on. Since there 
are two independent measurements, one each for the x and y 
directions, and each has the same statistical significance, the 
average of the two measurements is considered as the best es- 
timate ofR^"^^: 

R^^=^-(Rf^^+Rf^^) (23) 

The values of Gyis for each method and each scan are reported 
in TablefTOlin Appendix A. While the results of the second and 
third luminosity scans are compatible within statistical uncer- 
tainties, those of the first luminosity scan are lower by 2.7% to 
4.8% for all online algorithms, but are consistent for the offline 
track-based algorithms. These differences again indicate possi- 
ble systematic variations occurring between machine fills and 
are most likely to be caused by variations in the beam current 
cahbration (see Section l43T l. 

Figure|5](and TablefTOlin Appendix A) also report the spe- 
cific luminosity normalized to units of lO'^ protons per bunch 

X.pec = 1022(p/bunch)2^-^ (24) 

Because the emittance of Scan I was different from that of 
Scans II and III, the specific luminosity of that scan is not ex- 
pected to be the same as for the later scans. The agreement 
between algorithms within one scan is excellent. This agree- 
ment demonstrates that the variation in the measured value of 
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Gvis with scan number for a given algorithm is due to variations 
in the fitted value of R^^ rather than in the values obtained 
forZ. 

4.5 Systematic Uncertainties 

Systematic uncertainties affecting the luminosity and visible 
cross section measurements arise from the following effects. 

1. Beam intensities 

A systematic error in the measurement of the absolute bunch 
charge translates directly into an uncertainty on the lumi- 
nosity calibration. The accuracy of the bunch intensity mea- 
surement depends on that of the DCCT calibration. While 
laboratory measurements indicate an rms absolute scale un- 
certainty of better than 1.2%, the DCCT suffers from slow 
baseline drifts that are beam-, time- and temperature-depend- 
ent. These baseline offsets can only be determined with no 
beam in the LHC. 

For the fills under consideration, the DCCT baseline was 
measured before injection, and then again after dumping 
the beam. The DCCT-baseline determination is subject to 
magnetic and electronic drifts that translate into an rms un- 
certainty on the total circulating charge of ^ 1 . 15 x 10^ pro- 
tons. Conservatively combining the uncertainty on the ab- 
solute scale and on the baseline subtraction linearly yields 
a fractional uncertainty on the total charge ni(2) in beam 1 
(2) of 

C7(ni(2)) 1.15x10'^ 
^ ^ '^ = +0.012 (25) 

«1(2) «6"1(2) 

Treating the current-scale uncertainty as fully correlated 
between the two beams results in a total systematic error 
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Source 


Uncertainty on CTi.jj (% ) 


Beam Intensities 


10 


Length-Scale Calibration 


2 


Imperfect Beam Centring 


2 


Transverse Emittance Growth & Other Sources of Non-Reproducibility 


3 


fl Dependence 


2 


Fit Model 


1 


Total 


11 



Table 3. Summary of systematic uncertainties on the visible cross sections obtained from beam scans. Because CT,,/., is used to determine the 
absolute luminosity (see Equation[3](, these results are also the systematic uncertainty on the beam-scan based luminosity calibrations. 



of ±10% on the product of bunch currents for the running 
conditions summarized in Table 2. Because the baseline 
correction dominates the overall bunch-charge uncertainty, 
and because it drifts on the time scale of a few hours, these 
uncertainties are largely uncorrelated between the first (scan 
I) and the second (scans IIh-III) luminosity-calibration ses- 
sions. 

2. Length-Scale Calibration 

Fits to the beam size depend on knowledge of the relative 
displacement between the beams at each scan step. Thus, 
any miscaUbration of the beam separation length-scale will 
result in a mismeasurement of the luminosity. The desired 
nominal beam separation during beam scans determines the 
magnet settings of the closed orbit bumps that generate 
the beam separation. The only accelerator instrumentation 
available for calibrating the length-scale of the beam sepa- 
ration is the beam position monitor system. Unfortunately, 
the short-term stability and reliability of this system are 
not adequate to perform such a calibration. In contrast, the 



vertex resolution of the ATLAS Inner Detector provides 
a stable and precise method of calibration. These calibra- 
tions were done in dedicated scans where both beams were 
moved in the same direction first by +100 j^m and then 
by —100 fim from the nominal beam position, first in the 
horizontal and then in the vertical direction. The luminous 
beam centroid was determined using reconstructed primary 
vertices. In addition, the primary vertex event rate was mon- 
itored to ensure that the two beams remained centred with 
respect to each other. The calibration constants derived for 
the length-scale were (1.001 ±0.003) and (1.001 ±0.004) 
in the horizontal and vertical directions respectively, indi- 
cating that the scale associated with the magnet settings and 
that obtained from the ATLAS Inner Detector agree to bet- 
ter than 0.5%. The dominant source of uncertainty is the 
precision with which the two beams could be kept trans- 
versely aligned during the length-scale calibration scans. 
In addition, these scans consisted of only three points and 
extended to only ±100 jJ-m; therefore these data do not 



allow for studies of non-linearities, nor for checks of the 
calibration at the larger beam displacements used during 
the luminosity-calibration scans. Finally, if the transverse 
widths of the two beams happened to be significantly dif- 
ferent, the measured displacements of the luminous cen- 
troid at each scan point would not exactly reflect the av- 
erage displacement of the two beams. The combination of 
these effects results in an estimated systematic uncertainty 
of 2% on the length-scale calibration, in spite of the high 
precision of the calibration-scan data. 

3. Imperfect Beam Centring 

If the beams are slightly offset with respect to each other in 
the scan direction, there is no impact on the results of the 
luminosity scan. However, a deviation from zero separation 
in the transverse direction orthogonal to that of the scan 
reduces the rate observed for all the data points of that scan. 
The systematic uncertainty associated with imperfect beam 
centring has been estimated by considering the maximum 
deviation of the peak position (measured in terms of the 
nominal beam separation) from the nominal null separation 
that was calibrated through the re-alignment of the beams at 
the beginning of that scan. This deviation is translated into 
an expected decrease in rate and therefore in a systematic 
uncertainty affecting the measurement of the visible cross 
section. A systematic uncertainty of 2% is assigned. 

4. Transverse Emittance Growth and Other Sources of Non- 
reproducibiUty 

Wire-scanner measurements of the transverse emittances of 
the LHC beams were performed at regular intervals during 
the luminosity-scan sessions, yielding measured emittance 
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degradations of roughly 5-10% per beam and per plane be- 
tween the beginning and the end of the luminosity-calibration 
sessions l22l . This emittance growth causes a progressive 
increase of the transverse beam sizes (and therefore of Ex 
and Ey), leading to a ^ 2% degradation of the specific lu- 
minosity between the first and the last scan within one ses- 
sion. This luminosity degradation, in turn, should be re- 
flected in a variation over time of the specific rates 
and R^'^-^ (Eq.l23]). A first potential bias arises because if 
the time dependence of E^ and Ey during a scan is not taken 
into account, the emittance growth may effectively distort 
the luminosity-scan curve. Next, and because the horizon- 
tal and vertical scans were separated in time, uncorrected 
emittance growth may induce inconsistencies in computing 
the luminosity from accelerator parameters using Eq. |22] 
The emittance growth was estimated independently from 
the wire-scanner data, and by a technique that relies on the 
relationship, for Gaussian beams, between E, the single- 
beam sizes ai and 02 and the transverse luminous size Gl 
(which is measured using the spatial distribution of primary 
vertices) Il23ll : 

E = ^a^ + ai 

Here the emittance growth is taken from the measured evo- 
lution of the transverse luminous size during the fill. The 
variations in both E and r'^^^ (which should in princi- 
ple cancel each other when calculating the visible cross- 
section) were then predicted from the two emittance-growth 
estimates, and compared to the luminosity-scan results. While 
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the predicted variation of E between consecutive scans is 5 (CMS); the size of this correlated uncertainty remains under 
very small (0.3 - 0.8 /im) and well reproduced by the data, study, 
the time evolution of r'^^^ displays irregular deviations 
from the wire-scanner prediction of up to 3%, suggesting 

5 Internal Consistency of Luminosity 

that at least one additional source of non-reproducibility 

.... . Measurements 

is present. Altogether, these estimates suggest that a ±3% 
systematic uncertainty on the luminosity calibration be as- 

It is possible to test the consistency of the vdM calibrations by 

signed to emittance growth and unidentified causes of non- 
comparing the luminosities obtained using different luminosity 

reproducibiUty. 

detectors and/or algorithms. Figure |6] shows the instantaneous 

5. /i -Dependence of the Counting Rate 

luminosities obtained by various algorithms for Run 16288 

All measurements have been corrected for jj, dependent non- 
each normalized using the calibration extracted from its vdM 

linearities. Systematic uncertainties on the predicted count- 
scan data. The absolute luminosities agree to better than 2%; 

ing rate as a function of ji have been studied using Monte 

the relative luminosities track each other over time to within 

Carlo simulations, where the efficiency (or equivalently a,,n ) 

the statistical fluctuations. Over most of the 2010 pp run, LU- 

have been varied. For jj. <2 the uncertainty is estimated to 

CID_Event_OR was chosen as the preferred offline algorithm 

be < 2%, as illustrated in Fig.|2l 

because its pile-up correction was well-understood, its statis- 

6. Choice of Fit Model 

tical power was adequate and backgrounds for this algorithm 

For all methods, fits of the scan data to the default function 

were low. 

(double Gaussian with common mean plus constant back- 
Comparing the residual /i -dependence (if any) of the mea- 

ground for the online algorithms and double Gaussian for 

sured luminosity across multiple detectors and algorithms probes 

the background-free offline algorithms) have per degree 

the consistency of the pile-up correction procedures described 

of freedom values close to 1 .0, indicating that the fits are 

in Section 13.41 Figure |7] shows, for some of the LUCID and 

good. The systematic uncertainty due to this choice of fit 

MBTS algorithms, the raw counting rate as a function of the 

function has been estimated by refitting the offline data us- 

average number of inelastic interactions per BC measured by 

ing a cubic spline as an alternative model. The value of (7,,i.s 

LUCID.Event.OR using the prescription of Section U.4. II Non- 
changes by approximately 1%. 

linearities are apparent (as expected) for the LUCID_Event_AND, 
A summary of the systematic uncertainties is presented in Ta- LUCID.Event.OR and MBTS_Event_AND algorithms. If the para- 
ble[3] The overall uncertainty of 11% is dominated by the mea- metrizations of Section [34l are correct, however, then the ratio 



surement of the beam intensities. At least some portion of this ^ The bunch-by-bunch luminosity for LUCID_Event_OR averaged 
uncertainty is common to interactions points 1 (ATLAS) and over the full run is shown in Fig.[T] 
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of the luminosities determined using the different algorithms 
should be independent of /x. Figure [8] shows that the values of 
jj. obtained with the LUCID.EventJVND and MBTS.EventJVND 
algorithms remain within ±1% of that measured using the LU- 
CID_Event_OR algorithm over the range < ji < 2.5. Com- 
parisons of the LUCID.Event.OR and LUCID_Event_AND algo- 
rithms demonstrate agreement up to = 5, the highest value of 
jj. obtained during the 2010 LHC run. No results are presented 
beyond /i = 2.5 for the MBTS because during the correspond- 
ing data-taking period the short spacing between consecutive 
LHC bunches made the MBTS luminosity measurement unre- 
liable. Possible causes include the long duration of the analog 
pulse, saturation effects following large energy deposits, time 
jitter introduced by the electronics used at the time, and after- 
glow background. 

6 Comparison with Monte Carlo Generators 

Because the vdM method does not require knowledge of the 
inelastic cross section nor of the detector acceptance, the val- 
ues of (7,,;.s obtained from the beam scans can be used to test 
the accuracy of the predictions of Monte Carlo event genera- 
tors. Such predictions suffer from several theoretical uncertain- 
ties. First, because the pp inelastic cross section has not been 
measured at 7 TeV, the generators obtain (7,„g/ by extrapolating 
from lower energy. Results of this extrapolation depend on the 
functional form used. The PYTHIA and PHOJET generators, 
for example, predict values for (7,„e/ that differ by 6.6%. Sec- 
ond, the generators must separately model the non-diffractive 
(ND), single-diffractive (SD) and double-diffractive (DD) com- 
ponents of the cross section. There exists no unique prescrip- 
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Fig. 6. (a) ATLAS instantaneous luminosity for Run 162882, as mea- 
sured using several algorithms. Each curve is independently normal- 
ized using the vdM calibration obtained for that algorithm. The in- 
set at the bottom shows the ratio of the luminosity obtained with 
each algorithm to that obtained with LUCID.Event.OR. The statis- 
tical uncertainties for the online algorithms (LUCID.Event.OR, LU- 
CID.EventJ^ND and MBTS.EventJ^ND) are negligible. Statistical 
uncertainties for the offline algorithms (LAr.Tlming and ChPart) are 
displayed, (b) Comparison of the integrated luminosity obtained for 
Run 162882 for each of the algorithms shown above, together with the 
statistical uncertainties on the measurements. The dotted line shows 
the weighted mean of all the algorithms. The shaded band indicates a 
±2% deviation from that mean. 
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Fig. 7. Fraction of bunch crossings containing a detected event 
(N/Nfic) for several algorithms, as a function of ;ULuciD.Event.oR- 
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Fig. 8. Fractional deviation of the average value of jl obtained with the 
MBTS.EventJ\ND and LUCID_EventJ\ND algorithms with respect 
to the LUCID_Event_OR algorithm as a function of fl obtained with 
LUCID_EvenLOR. 

tion for classifying events as diffractive or non-diffractive and 
no calculation of the cross sections from first principles. Typ- 
ical uncertainties associated with such classifications are il- 
lustrated in Table |4] The fraction of (7,„p/ corresponding to 
ND events is 68% in PYTHIA and 81% in PHOJET, while 
the DD fractions are 13% and 5% respectively. Finally, there 
are significant uncertainties on the modeling of the predicted 



multiplicity-, pj- and Tj- distributions for particles produced 
in soft pp interactions, particularly for the poorly constrained 
diffractive components. Differences in these distributions will 
affect the efficiencies for events to pass the selection criteria of 
a specific luminosity algorithm. 

Within the framework of Monte Carlo generators, CJv,i is 
calculated using the expression 



<7,.,s — £ND(yND + £SD<^SD + ^DDdDD 



(27) 



where Eprocess are the efficiencies and Gprocess the cross sections 
for the individual inelastic processes (ND, SD and DD). Table|5] 
shows the predicted efficiencies for observing ND, SD and DD 
events using either PYTHIA (with the default ATLAS MC09 
tune ll24l ) or PHOJET, for some of the algorithms described in 
Section [3] In general, the PHOJET predictions are about 15% 
to 20% higher than those obtained with PYTHIA. One excep- 
tion is LUCID_Event_AND which is less sensitive to the diffrac- 
tive processes: here the two generators agree to within 5% over- 
all. Additional systematic uncertainties on these predictions, 
associated with the modeling of the detector response in the 
simulation, are algorithm- and trigger-dependent and vary from 
2.2% for MBTS.Event.OR to 6% for 
LUCID.EventJ^ND. 

As noted in Section 14.41 there is a systematic difference 
between the values of Ovis obtained from the first scan and 
those based on the second and third scans. In reporting our 
best estimate of the measured visible cross sections, we chose 
to average the results of the first scan with the average of the 
second and third scans. Comparisons of the vdM scan mea- 
surements with the Monte Carlo predictions are presented in 
Tableland Fig. |9l For a given event generator, the compar- 



Cross Section &t^/s = l TeV 


Process 


PYTHIA 


PHOJET 




(mb) 


(mb) 


non-dii&active (ND) 


48.5 


61.6 


single-diffractive (SD) 


13.7 


10.7 


double-diffractive (DD) 


9.3 


3.9 


Total: 


71.5 


76.2 



Table 4. Predicted inelastic pp cross sections dX yfs = 1 TeV for 
PYTHIA and for PHOJET. A small (~ 1 mb) contribution from 
double-pomeron processes ("central diffraction") was not included in 
the PHOJET cross section. 

isons exhibit an RMS spread of 4 to 5%; on the average, the 
PYTHIA (PHOJET) predictions are 15% (33%) higher than 
the data. Given the 1 1% systematic uncertainty on the vdM cal- 
ibration, which is correlated across all algorithms, PYTHIA 
agrees with the data at the level of 1.5 a, while PHOJET and 
the data deviate at the 3(7 level. 

7 Conclusions 

Measurements of the LHC luminosity have been performed by 
ATLAS in proton-proton collisions at ^/S = 7 TeV using mul- 
tiple detectors and algorithms. The absolute luminosity cali- 
brations obtained using beam- separation scans suffer from a 
±11% systematic uncertainty, that is dominated by the uncer- 
tainty in the bunch intensities and is therefore highly correlated 
across all methods. For a given bunch luminosity, i.e. for a fixed 
value of jU (the average number of inelastic pp interactions per 
crossing), the absolute luminosities obtained using different de- 
tectors and algorithms agree to within ±2%. In addition, the 
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luminosities from these methods track each other within better 
than 2% over the range < /i < 2.5. The visible cross sections 
obtained from the beam scan cahbrations also have a system- 
atic uncertainty of 11% and are lower than those predicted by 
PYTHIA (PHOJET) by about 15 % (33%). 
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LUCID.Event.OR 


LUCID.EventJVND 


Process 


Efficiency (% ) 


Efficiency (% ) 




PYTHIA MC09 


PHOJET 


PYTHIA MC09 


PHOJET 


ND 


79.7 


73.7 


30.8 


24.9 


SD 


28.7 


44.3 


1.3 


2.4 


DD 


39.9 


62.0 


4.3 


14.6 


Ovis (mb) 


46.4 


53.1 


16.0 


17.0 





MBTS.Timing 


LAr.Timing 


Process 


Efficiency (% ) 


Efficiency (% ) 




PYTHIA MC09 


PHOJET 


PYTHIA MC09 


PHOJET 


ND 


97.4 


97.9 


96.0 


94.3 


SD 


41.3 


44.3 


21.4 


27.9 


DD 


50.8 


68.1 


25.9 


53.6 


(J,,is (mb) 


57.6 


67.8 


51.9 


63.2 





ChPart 


PrimVtx 


Process 


Efficiency (% ) 


Efficiency (% ) 




PYTHIA MC09 


PHOJET 


PYTHIA MC09 


PHOJET 


ND 


85 


80 


97.8 


99.2 


SD 


36 


36 


43.9 


56.9 


DD 


36 


41 


47.8 


70.7 


CT,,,-.5 (mb) 


45.7 


54.7 


57.9 


70.0 



Table 5. Efficiencies at = 7 TeV for several of the luminosity methods described in Section[3] The predicted visible cross sections (7,.,j are 
obtained using Equation|27] the efficiencies in the present table and the cross sections in Tablej?) 
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Algorithm 


(mb) 


_PYTHIA 

vis 

(mb) 


^PYTHIA 


JJ.PHOJET 

(mb) 


jyPHOJET 






LUCID.EventJ\ND 


12.4±0.1 


16.0±0.8 


1.29±0.07 


17.0±0.9 


1.37 ±0.07 


LUCID.EvenLOR 


40.2 ±0.1 


46.4 ±2.8 


1.15±0.07 


53.1±3.2 


1.32 ±0.08 


MBTS.EventJ\ND 


51.9±0.2 


58.4±1.5 


1.13±0.03 


68.7± 1.8 


1.32±0.03 


MBTS.EvenLOR 


58.7±0.2 


66.6 ±1.5 


1.13±0.03 


73.7±1.6 


1.26 ±0.03 


MBTS.Timing 


50.4 ±0.2 


57.6±1.3 


1.14±0.03 


67.8±1.8 


1.35 ±0.04 


PrimVtx 


53.6±0.2 


57.9±1.3 


1.08 ±0.03 


70.0±1.6 


1.31±0.03 


ChPart 


42.7 ±0.2 


45.7 ±1.7 


1.07 ±0.04 


54.7±2.0 


1.28 ±0.05 


LAr.Timing 


46.6 ±0.2 


51.9±2.3 


1.11±0.05 


63.2±2.9 


1.36±0.06 



Table 6. Comparison of the visible cross sections determined from beam scans {a^is) to the predictions of the PYTHIA and PHOJET Monte 
Carlo generators. The ratio of prediction to measurement is also shown. The errors affecting the measured visible cross sections are statistical 
only. The errors on the PYTHIA and PHOJET visible cross sections are obtained from the systematic uncertainty associated with modeUng 
the detector response. These uncertainties are fully correlated, row by row, between PYTHIA and PHOJET; they are fully correlated between 
the two LUCID algorithms, and highly correlated for the five MBTS-triggered algorithms (MBTS J\ND, MBTS.OR, MBTS.timing.Event, 
PrimVtx_Event and Ch Part-Event). The fully correlated 11% systematic uncertainty on visible cross sections, that arises from the vdM 
calibration, is not included in the errors listed in this table. 

LAS Tier-1 faciUties at TRIUMF (Canada), NDGF (Denmark, 
Norway, Sweden), CC-IN2P3 (France), KIT/GridKA (Germany), 
INFN-CNAF (Italy), NL-Tl (Netherlands), PIC (Spain), ASGC 
(Taiwan), RAL (UK) and BNL (USA) and in the Tier-2 facih- 
ties worldwide. 



30 



ATLAS 



ATLAS 



, -> 
1 ■a 



• LUCID_Event_AND 

■ LUCID_Event_OR 

O MBTS_AND 

□ MBTS_OR 

A MBTS_timing_Event 

PrimVtx Event 

^ ChPart_Event 

LAr timing Event 
in ZDC_A 
X ZDC_C 



• LUCID_Event_AND 
■ LUCID_Event_OR 
O MBTS_AND 

□ MBTS_OR 

A MBTS timing Event 

PrimVtx Event 

* ChPart_Event 
LAr timing Event 

yn ZDC_A 
X ZDC_C 



57 



58 59 60 



(a) 



61 



62 



60 



61 



62 



63 



(b) 



64 



65 66 

£y [^im] 



6.5 



7.5 



• LUCID_Event_AND 
■ LUCID_Event_OR 
O IVIBTS_AND 

□ IUIBTS_OR 

A IvlBTS_timing_Event 

PrimVtx Event 

* CtiPart_Event 
LAr timing Event 

^ ZDC_A 

X ZDC_C 



8.5 9 9.5 



ATLAS 

• LUCID_Event_AND 
■ LUCID_Event_OR 
O IVIBTS_AND 

□ IVIBTS_OR 

A IvlBTSJiming Event 

^ PrimVtx Event 

* ChPart_Event 

^ LAr timing Event 



ZDC A 



1 1.5 



2.5 



3.5 



4.5 5 



(c) (d) 
Fig. 4. Fit results for the values of (a) Lx, (b) Ey, (c) xq and (d) yo obtained using different luminosity algorithms during Scan II. The dashed 
vertical line shows the unweighted average of all the algorithms. The shaded bands indicate ±0.5% deviations from the mean for (a) and (b) 
and ±0.1 /im deviations from the mean for (c) and (d). In all cases, the uncertainties on the points are the statistical errors reported by the vdM 
fit. Uncertainties for different algorithms using the same detector are correlated. 



31 



ATLAS 



• LUCID_Event_AND 
■ LUCID_Event_OR 
O MBTS_AND 
□ MBTS_OR 
A MBTS_timing_Event 
PrimVtx Event 
ChPart_Event 
^ LAr timing Event 
in ZDC_A 
X ZDC_C 



ATLAS 

• LUCID_Event_AND 
■ LUCID_Event_OR 
O MBTS_AND 
□ MBTS_OR 
A MBTS_timing_Event 
PrimVtx Event 
ChPart_Event 
^ LAr timing Event 
in ZDC_A 
X ZDC_C 



4.8 4.85 4.9 4.95 5 5.05 5.1 4.8 4.85 4.9 4.95 5 5.05 5.1 

Lspec [ 1 cm-2s-i bunch"^ (10" p)-^ ] 1^^^^ [ 1 0^^ cm-^s-^ bunch"^ (10" p)-^ ] 

(a) (b) 

Fig. 5. Comparison of the specific luminosities obtained using various luminosity algorithms for (a) Scan II and (b) Scan III. The dashed lines 
show the unweighted average of all algorithms; the shaded band indicates a ±0.5% variation from that mean. The uncertainties on the points 
are the statistical errors reported by the vdM fit. Uncertainties for different algorithms using the same detector are correlated. 
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(a) (b) 
Fig. 9. Comparison of the measured values of CT,,,^ for several algorithms to the (a) PYTHIA and (b) PHOJET predictions. The errors on 
the points are the systematic uncertainties due to possible inaccuracies in modeling the detector response. The uncertainties for different 
algorithms using the same detector are correlated. The 11% uncertainty on the vdM calibration of the luminosity, which is 100% correlated 
among algorithms, is not included in the error bars. 
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A Fits to Beam Scan Data 

This appendix presents results of the fits to vdM scan data for all scans and all algorithms. 



Algorithm 


Mean Position 


E 


Background 




X^/DOF 






(jum) 


(Hz) 


(Hz) 




Horizontal Scan I 


LUCID_EventJ\ND 


-1.12± 0.46 


47.40 ± 0.56 


0.01 ± 0.04 


75.6 ± 1.1 


0.9 


LUCID.Event.OR 


-1.58± 0.25 


47.27 ± 0.29 


0.06 ±0.04 


247.8 ± 2.0 


0.5 


MBTSJ^ND 


-1.85± 0.25 


47.33 ± 0.25 


0.03 ± 0.04 


319.0 ± 2.3 


0.8 


MBTS.OR 


-2.05± 0.24 


47.30 ± 0.26 


1.01 ±0.11 


361.7 ± 2.6 


1.0 


MBTS.timing.Event 


-1.66 ±0.26 


47.05 ±0.26 


N/A 


306.8 ±1.6 


1.0 


PrimVtx.Event 


-1.7± 0.2 


47.26 ±0.25 


N/A 


329.7 ±1.6 


0.8 


ChPart-Event 


-1.67±0.3 


47.3 ±0.3 


N/A 


253.2 ±1.6 


0.8 


LAr_timing_Event 


-1.44 ±0.27 


47.0±0.3 


N/A 


290.6±1. 9 


0.5 


BCM.Event_OR 


-2.33 ±1.42 


47.27 (fixed) 


7.5 ±0.20 


26.98 ±0.89 


0.6 


Vertical Scan I 


LUCID.EventJ\ND 


-5.04± 0.50 


55.52 ± 0.59 


0.05 ± 0.03 


75.8 ± 1.0 


0.8 


LUCID.Event.OR 


-5.23± 0.28 


55.28 ± 0.33 


0.16 ±0.06 


246.2 ± 1.9 


1.1 


MBTSJiiND 


-5.24± 0.28 


55.73 ± 0.30 


0.10 ±0.06 


318.5 ± 2.3 


1.2 


MBTS.OR 


-5.25± 0.26 


55.82 ± 0.28 


1.08 ±0.12 


359.2 ± 2.5 


1.2 


MBTS.timing_Event 


-5.53 ±0.30 


56.32 ±0.29 


N/A 


297.8 ±1.4 


2.1 


PrimVtx_Event 


-5.17 ±0.26 


56.28 ±0.30 


N/A 


323.0 ±1.5 


1.1 


ChPart.Event 


-5.61± 0.35 


56. 1 ± 0.4 


N/A 


249. 3 ± 1.6 


1.4 


LAr_timing_Event 


-5.11 ±0.31 


56.2 ±0.4 


N/A 


280.6±1. 8 


2.1 


BCM.Event_OR 


-3.63±1.51 


55.28 (fixed) 


7.5 ±0.20 


27.3 ±0.8 


0.7 



Table 7. Summary of the relevant fit parameters for the Beam Scan 1. For offline algorithms, the rates have been corrected for trigger prescales. 
Because the rates in the BCM were low, the value of E used for the BCM was fixed to that obtained from the LUCID_Event_OR. No results 
are presented for the ZDC, since the constant fraction discriminators used for the ZDC measurements were installed later in the run. 
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Algorithm 


Mean Position 


E 


Background 




X^/DOF 




(jum) 


{m) 


(Hz) 


(Hz) 




Horizontal Scan n 


LUCID.EventJ!iND 


7.65 ± 0.25 


58.78 ± 0.16 


-0.02± 0.06 


265.4 ± 3.0 


1.8 


LUCID.Event.OR 


7.41 ± 0.14 


58.76 ± 0.08 


0.07 ±0.12 


858.9 ± 2.5 


2.0 


MBTSJVND 


7.28 ± 0.13 


59.06 ± 0.09 


-0.28± 0.16 


1107.3 ±3.1 


0.9 


MBTS_OR 


7.30 ±0.13 


58.93 ± 0.09 


1.04 ± 0.25 


1253.1 ±3.6 


1.2 


MBTS_timing_Event 


7.44 ±0.22 


58.71 ±0.23 


N/A 


1087.0 ±4.1 


1.3 


PrimVtx_Event 


7.56 ±0.20 


58.63 ±0.21 


N/A 


1133.0±4.0 


1.1 


ChPart_Event 


7.42 ±0.34 


58.5 ±0.2 


N/A 


869.1 ±4.2 


1.1 


LAr_timing_Event 


7.41 ±0.28 


58.2±0.3 


N/A 


997.5 ±5.6 


1.6 


BCM_Event_OR 


6.54 ±0.59 


58.76 (fixed) 


0.313 ±0.083 


89.00 ±0.95 


0.9 


ZDCJV 


6.98 ±0.22 


59.05 ±0.12 


0.09 ±0.14 


380.7 ±1.8 


1.1 


ZDC_C 


6.88 ±0.24 


58.74 ±0.1 9 


0.32 ±0.10 


370.57 ±2.0 


0.8 


Vertical Scan n 


LUCID.EventJ^ND 


1.99 ±0.27 


62.75 ±0.19 


-0.21± 0.14 


253.8 ± 2.9 


1.6 


LUCID.Event.OR 


1.99 ±0.16 


62.37 ± 0.16 


0.13 ±0.13 


825.3 ±3.1 


0.8 


MBTSJVND 


2.17 ±0.15 


62.18 ±0.16 


0.30 ±0.15 


1068.9 ± 3.9 


0.9 


MBTS.OR 


2.11 ±0.15 


62.13 ± 0.15 


1.70 ±0.20 


1207.6 ± 4.2 


1.0 


MBTS_timing_Event 


2.22 ±0.24 


62.61 ±0.27 


N/A 


1038.0 ±3.8 


1.5 


PrimVtx.Event 


2.09 ±0.21 


62.48 ±0.25 


N/A 


1081.0±3.6 


0.9 


ChPart_Event 


2.27± 0.36 


62.3 ±0.3 


N/A 


841. 2 ±4.1 


1.1 


LAr_timing_Event 


2.50 ±0.29 


62.7 ±0.4 


N/A 


950.6 ±5.4 


3.0 


BCM_Event_OR 


1.85 ±0.63 


62.37 (fixed) 


0.429 ±0.079 


85.53 ±0.89 


1.2 


ZDCJ\ 


2.54 ±0.25 


62.00 ±0.27 


0.45 ±0.12 


368.9 ±2.3 


1.1 


ZDC.C 


2. 15 ±0.25 


62.38 ±0.28 


0.34±0.12 


355.9±2.3 


0.8 



Table 8. Summary of the relevant fit parameters for the Beam Scan II. For offline algorithms, the rates have been corrected for trigger prescales. 
Because the rates in the BCM were low, the value of X used for the BCM was fixed to that obtained from the LUCID.Event.OR. 
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Algorithm 


Mean Position 


E 


Background 




X^/DOF 






{m) 


(Hz) 


(Hz) 




Horizontal Scan in 


LUCID.EventJVND 


5.48 ± 0.26 


58.94 ±0.19 


0.04 ±0.13 


266.8 ± 3.0 


1.2 


LUCID.Event.OR 


5.66 ± 0.15 


58.57 ± 0.18 


0.42 ±0.10 


856.8 ± 3.3 


2.1 


MBTSJVND 


5.59 ± 0.14 


58.88 ± 0.10 


0.15 ±0.14 


1102.5 ±3.2 


2.3 


MBTS_OR 


5.59 ± 0.14 


58.87 ±0.10 


1.20 ± 0.30 


1244.4 ± 3.9 


2.5 


MBTS_timing_Event 


6.02 ±0.22 


59.05 ±0.23 


N/A 


1074.0 ±4.0 


0.95 


PrimVtx_Event 


5.95 ±0.20 


59. 14 ±0.23 


N/A 


1120.0±3.8 


1.4 


ChPart_Event 


6.03 ±0.33 


59.3 ±0.2 


N/A 


869.6 ±4.2 


1.1 


LAr_timing_Event 


6.15±0.28 


59.1 ±0.3 


N/A 


981.7±6.6 


1.4 


BCM_Event.OR 


6.36 ±0.60 


58.57 (fixed) 


0.23 ±0.11 


89±1 


1.25 


ZDCJV 


5.38 ±0.22 


59.15 ±0.36 


0.28 ±0.18 


373.6±3.1 


1.3 


ZDC_C 


5.67 ±0.23 


59.01 ±0.15 


0.13 ±0.10 


366. 7 ± 1.8 


1.7 


Vertical Scan m 


LUCID.EventJVND 


-0.01±0.27 


62.21 ± 0.30 


-0.03± 0.08 


259.9 ± 2.9 


0.9 


LUCID.Event.OR 


0.08 ±0.16 


62.06 ± 0.16 


0.23 ±0.12 


830.2 ±3.1 


0.8 


MBTSJVND 


0.04 ± 0.15 


62.09 ± 0.16 


0.15 ±0.15 


1075.6 ± 3.9 


1.2 


MBTS.OR 


0.06 ± 0.15 


62.09 ± 0.15 


1.65 ± 0.22 


1214.5 ± 4.2 


1.1 


MBTSJming.Event 


-0.16 ±0.24 


61.45 ±0.30 


N/A 


1056.0±4.0 


1.4 


PrimVtx.Event 


-0.06 ±0.21 


61.83±0.27 


N/A 


1102.0±3.7 


1.4 


ChPart_Event 


-0.32± 0.36 


61.5±0.3 


N/A 


840.6 ±4.1 


0.9 


LAr_timing_Event 


-0.53 ±0.30 


61.7±0.5 


N/A 


951. 1±6 .2 


3.6 


BCM_Event.OR 


0.3 ±0.64 


62.06 (fixed) 


0.17 ±0.08 


86.2 ±1 


1.56 


ZDCJ\ 


-0.04 ±0.25 


62.36 ±0.28 


0.17 ±0.10 


367.9 ±2.3 


1.2 


ZDC.C 


-0.03 ±0.25 


62.26 ±0.30 


0.42±0.10 


358.3±2.3 


0.8 



Table 9. Summary of the relevant fit parameters for the Beam Scan III. For offline algorithms, the rates have been corrected for trigger prescales. 
Because the rates in the BCM were low, the value of X used for the BCM was fixed to that obtained from the LUCID.Event.OR. 
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Method 


Scan Number 










mb 


(lO^^cm-^s-i) 


LUCID_EventJ\ND 


1 


12.15±0.14 


6.80 ±0.08 




2 


12.55±0.10 


4.85 ±0.03 




3 


12.73±0.10 


4.88 ±0.03 


LUCID.EvenLOR 


1 


39.63 ±0.32 


6.85 ±0.06 




2 


40.70±0.13 


4.88±0.01 




3 


40.77 ±0.14 


4.92 ±0.02 


MBTS.EventJ!iND 


1 


51.14±0.39 


6.78 ±0.05 




2 


52.59±0.16 


4.87 ±0.01 




3 


52.64±0.16 


4.90 ±0.02 


MBTS.Event.OR 


1 


57.83±0.43 


6.79 ±0.05 




2 


59.47 ±0.18 


4.89 ±0.01 




3 


59.43 ±0.25 


4.90 ±0.02 


MBTS.Timing 


1 


49.28 ±0.31 


6.76 ±0.05 




2 


5 1.64 ±0.23 


4.87 ±0.03 




3 


5 1.29 ±0.24 


4.93 ±0.03 


PrimVtx 


1 


53.48 ±0.29 


6.73 ±0.05 




2 


53.64 ±0.22 


4.89 ±0.03 




3 


53.78 ±0.23 


4.89 ±0.02 


ChPart 


1 


42.61 ±0.26 


6.74 ±0.05 




2 


42.84±0.21 


4.91 ±0.02 




3 


42.93 ±0.21 


4.91 ±0.03 


LAr.Timing 


1 


46.43 ±0.31 


6.78 ±0.05 




2 


46.98 ±0.27 


4.91 ±0.02 




3 


46.63 ±0.30 


4.91 ±0.03 


ZDCJi 


2 


18.12±0.09 


4.89 ±0.02 




3 


17.78±0.13 


4.85 ±0.04 


ZDC.C 


2 


17.56±0.09 


4.88 ±0.02 




3 


17.39±0.11 


4.86 ±0.03 



Table 10. Measurements of the visible cross section and peak specific luminosity for all algorithms that have been calibrated using the vdM 
scan data for each of the three beam scans. The uncertainties reported here are statistical only. The emittance during Scan I was different from 
that during Scans 11 and III, so the specific luminosity in that first scan is not expected to be the same. No results for Scan I are presented for 
the ZDC, since the constant fraction discriminators used for the ZDC measurements were installed later in the run. Because the rates in the 
BCM were low, the value of E used for the BCM was fixed to that obtained from the LUCID.Event.OR, so no measurement of the specific 
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N. Ghodbane^^, B. Giacobbel'^^ S. Giagu^^^a.^ib^ y Giakoumopoulou^ V. Giangiobbei22a,i22b^ p Qianotti^'', B. Gibbard^^, 
A. Gibsonl5^ S.M. Gibson^^, G.F. Gieraltowski^, L.M. Gilbert^i^ M. Gilchriesei'*, O. Gildemeister^^, V. Gilewsky^i, 
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O. Gutzwiller"^ C. Guyot^^^ C. Gwenlan"^ C.B. Gwilliam^^^ A. Haas^^s, S. Haas^^, C. Haber^^, R. Hackenburg^^, 

H.K. Hadavand^^, D.R. Hadley", P Haefner^^, R. Hartel^^, R Hahn^^ S. Haider^^ Z. Hajduk^*, H. Hakobyani"^^, J. Haller^'^, 

K. Hamacher"^^ ^ Hamilton'*^, S. Hamiltoni^i, H. Han32'>, L. Han32b, K. Hanagaki"^, M. Hancei^", C. Handel^\ P Hanke^*^, 
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R. Hertenberger^^ L. Hervas^^, N.P. Hesseyi''^ A. ffidvegil'*^^ E. ffigon-Rodriguezi^"^, D. ffiU^-*, J.C. ffitf^, N. ffiU^, 
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